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ABSTRACT 

The global structure of optically thin advection dominated accretion flows 
which are composed of two-temperature plasma around black holes is calculated. 
We adopt the full set of basic equations including the advective energy 
transport in the energy equation for the electrons. The spectra emitted by the 
optically thin accretion flows are also investigated. The radiation mechanisms 
which are taken into accout are bremsstrahlung, synchrotron emission, and 
Comptonization. The calculation of the spectra and that of the structure of the 
accretion flows are made to be completely consistent by calculating the radiative 
cooling rate at each radius. As a result of the advection domination for the ions, 
the heat transport from the ions to the electrons becomes practically zero and 
the radiative cooling balances with the advective heating in the energy equation 
of the electrons. Following up on the successful work of Narayan et al. (1995), 
we applied our model to the spectrum of Sgr A*. We flnd that the spectrum of 
Sgr A* is explained by the optically thin advection dominated accretion flow 
around a black hole of the mass Mbh = 10^ M©. The parameter dependence of 
the spectrum and the structure of the accretion flows is also discussed. 

Subject headings: accretion, accretion disks — black hole physics — radiation 
mechanisms: non-thermal — Galaxy: center 
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1. Introduction 

Ever since the pioneering studies of steady thin accretion disks by Shakura & Sunyaev 
(1973, hereafter SS), the model of thin accretion disks has been apphed successfully to low 
energy emission from astrophysical objects, such as dwarf novie and pre-main-sequence 
stars. However, the model has been less successful in modeling of high energy emission 
from Galactic black hole candidates and active galactic nuclei which are considered to be 
powered by accreting black holes. The main problem lies in the fact that the thin accretion 
disks can not reproduce the observed spectra of such systems. The thin accretion disk 
model which SS originally proposed assumes that the disk is optically thick in the vertical 
direction and radiates the energy generated by the viscosity locally. This model predicts 
that the generated spectrum is multi-colored black body, which cannot explain the observed 
power-law spectra of X-rays from AGNs and Galactic black hole candidates even though it 
can explain the UV bump of the AGN or the soft state spectrum of the Galactic black hole 
candidates. More fundamentally, although it is generally believed that QSOs and Seyfert 
galaxies are powered by the gas accretion onto a super-massive black hole of the mass 
Mbh ^ IO^Mq, the thin accretion disks onto such super-massive black holes are too cool to 
generate high energy photons which are observed in many QSOs and Sayfert galaxies. 

The model which is investigated by Shapiro, Lightman & Eardley (1976, hereafter 
SLE) is quite attractive in that the accreted gas is optically thin and is much hotter 
than that in the SS solution and is hot enough to produce high energy photons. SLE 
considered two-tcmpcraturc plasma with ions being much hotter than electrons, which 
enabled quantitative studies of the non-blackbody spectra. SLE model has been applied to 
explain the spectrum of X-ray binaries and active galactic nuclei successfully. However, it 
is known that the optically thin hot accretion disk which is considered in SLE is thermally 
unstable (Piran 1978). If the accretion disk is heated up, then the disk expands and the 
density decreases, so that the bremsstrahlung cooling rate decreases. The reduced cooling 
then causes the gas to become even hotter, leading to a runaway thermal instability. For 
this reason it is not likely that such hot accretion disks exist in real systems for much longer 
than the thermal timescale. 

The models introduced so far are local solutions in the sense that the heat generated via 
viscosity is locally radiated away efficiently, which corresponds to neglecting the advective 
energy transport in the energy equations. When plasma cannot emit radiation efficiently, 
the heat generated via viscosity is advected inwards as the internal energy of the plasma. 
Abramowicz ct al. (1988) investigated the effect of advection term in their "slim disk" 
model in the optically thick case, and made it clear that there exists advection dominated 
branch where the viscous heating is balanced with the advection term rather than the 
radiative cooling term at high mass accretion rates. 

The optically thin advection dominated solution at low mass accretion rates is studied 
intensively by Abramowicz et al. (1995) and Narayan &Yi (1995) (see also Ichimaru 1977, 
Matsumoto et al. 1985). Although they claimed that optically thin advection dominated 
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solution is thermally stable for the long wavelength perturbations, Kato et al. (1996) 

showed the possibilities of the instability against the short wavelength perturbations. 
Manmoto et al. (1996) demonstrated that such an instability is favored for explaining 
rapid X-ray fluctuations from Galactic black hole candidates and does not affect the global 
stabihty of the accretion flows. 

The optically thin advection dominated accretion flows are applied to explain the 
observed spectra of accreting black holes. Narayan & Yi (1995) investigated the self-similar 
solution which is used later to calculate the spectrum from several low-luminosity accreting 
systems with great success. As a next step, the calculations of full global steady solutions 
were awaited for further investigations. Chen et al. (1997) solved optically thin advection 
dominated solution globally, but their solutions are those of one-temperature plasma and 
do not include detailed radiation mechanisms. Narayan et al. (1997) derived global solution 
for optically thin advection dominated accretion flows and showed that the self-similar 
solution is a good approximation at the radius far enough from the outer and the inner 
boundaries. This means that the spectra which are derived by using the self-similar 
solution may be modified, because considerable amount of photons may come from the 
hot region near the inner boundary. Narayan et al. (1997) did calculate the spectra from 
two-temperature accretion flows with their global solutions, but they treated the electron 
energy equations locally, neglecting the effect of the electron advection. Nakamura et al. 
(1996) were the first to solve the energy equations for ions and electrons and obtained the 
global two-temperature advection dominated solutions, and showed that the temperature 
profiles which arc crucial to the generated spectra arc largely modified when the effect of 
electron advection is taken into account. However, Nakamura et al. (1996) focused their 
attentions on the structure of the optically thin accretion disks and did not investigate the 
spectra from the disks. 

The study on the spectrum from optically thin advection dominated accretion flows 
with full global treatment of the basic equations is yet to be done. Thus we are motivated to 
consistently solve full set of equations including the energy equation for ions and electrons 
with detailed radiation mechanisms and obtain the spectra from the optically thin accretion 
flows. 

In section 2, we present the physical assumptions and the basic equations of our model. 
We show the results of our calculations in section 3. We then apply our model to Sgr A* 
(the central core of our Galaxy) in section 4. We conclude in section 5 with a summary and 
discussion. 



2. Accretion Flow Model 
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2.1. Physical Assumptions 

We consider an optically thin, steady axisymmetric accretion flow around a black 
hole. To investigate the spectra generated by the optically thin gas flows, we discuss gas 
dynamics in the context of two-temperature plasma. Assuming the existence of randomly 
oriented magnetic flelds which possibly originate from the turbulence in the gas flow, we 
take total pressure p to be 

P = Pgas+Pmag, (1) 

where Pgas is the gas pressure, Pmag is the magnetic pressure. We neglect the radiation 
pressure in this paper because the optically thin accretion flows we consider is always gas 
pressure dominated. We take the ratio of gas pressure to the total pressure as a global 
parameter which we designate as (3. Technically we need to solve the magnetic fleld 
self-consistently with the gas dynamics, but it is beyond the scope of this paper to treat full 
magneto-hydrodynamics equations. Due to the two-temperature assumption, we write Pgas 
as k k 

Pgas = f3p = Pi + Pe = — Ti + — Te. (2) 

Hi mn /ie mn 

Here and hereafter subscripts i and e indicate the quantities for ions and electrons, 
respectively. In eq. (|^), T is the temperature, p is the density, and fi is the mean molecular 
weight which is given by 

fii = 1.23, /ie = 1.14, (3) 

where numerical values correspond to the cosmic abundance. We estimate magnetic fleld B 
via magnetic pressure by the following equation: 

Pmag = (1 -/?)P= — ■ (4) 



2.2. Basic Equations 

We adopt cylindrical coordinate system (r,Lp,z) to describe axisymmetric ( ^ = ) 
accretion flows. Basic equations which describe the dynamics of the accretion disks are 
those of mass conservation, Euler equations which comprises three spatial components, and 
the energy equations. Mass conservation gives 

d d 

— {rpvr) + r— {pv^) = 0, (5) 

where Vr and are the radial and the vertical velocity respectively. The radial component 
of Euler equation is 

d f 2\ ^ / ( dil) 2\ 

— [rpv,) + r— {pvrv,) = -p\r—-vJ- r— , (6) 



where v^p is the azimuthal velocity and is the potential energy. To simulate general 
relativistic effects, we adopt pseudo-Newtonian potential 



= -GMbh/{R - rg) (7) 

with Mbh being the mass of the black hole, the Schwarzschild radius, and R = [r"^ + z^Y/"^ 
the distance from the central black hole. This potential is known to represent the dynamical 
aspects of general relativistic effects quite well for r > 2rg (Paczynski & Wiita 1980) and 
greatly simplifies the basic equations. Technically we need to solve the basic equations in 
the Schwarzschild or Kerr metric especially for the photon propagation. However solving 
the equations including the photon propagation in the full relativisitc metric is out of the 
scope of current paper. Azimuthal component of Euler equation is the conservation of 
angular momentum which gives 

^ {r'^fyvrv^) + {rfw^v,) = ^ (rV,^) . (8) 

Here r^^ is the ryj-component of the stress tensor. According to the conventional 
aZp/ia-prescription of sheer viscosity, r^^ can be written as 

where is the angular velocity and Ok is the Keplcrian angular velocity on the equatorial 
plane. In our paper, we take r^^p to be simply proportional to the local pressure p: 

= -ap, (10) 

where a is the dimensionless viscosity parameter which in general is considered to be around 
0.01 ~ 0.1. Two-tcmpcrature assumption requires two energy equations i.e. one for ions 
and one for electrons, both of which contain advective energy transport terms. The energy 
equation for each species is 

^ / dsi dsi\ d^l ^ . . 

P^i y'r-g:^ + ^--^ J = - Aie, (11) 

(Os Os \ 

Vr-g^ + V,-^ j=Xie- q;^, (12) 

where s is the specific entropy, A^e is the volume energy transfer rate from ions to electrons, 
and is the volume radiative cooling rate. We assume that the energy is transferred from 
ions to electrons via Coulomb collisions only. Stepney & Guilbert (1983) give an explicit 
expression: 

\ _ 1 OK 3me^ ^ „ „ (kT,-kT,,:) 1 A 
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where the K^s are modified Bessel functions, In A = 20 is tlie Coulomb logaritfim, and 
6 = kT/nieC^ is tlie dimensionless temperature. Note tliat tfie energy equations assume 
tliat tlie ions and the electrons are in thermal equilibrium by some mechanisms. The 
thermodynamic relations are 

pTidsi = — ^ ( dpi - i—dp \ , (14) 

pTedSe = ( dpe - 'J^dp ) , (15) 

where 7 = 5/3 is the adiabatic index. In the energy equations we have assumed viscous 
heating acts only on ions because an ion particle is much heavier than an electron particle. 

Now we have full set of basic equations which describe the steady accretion flow around 
a black hole. However, it is rather a difficult problem to solve above partial differential 
equations with respect to r and z plus radiative transfer equations to obtain global steady 
solutions. To lessen the number of variables, we fix the vertical structure before we solve 
the full spatial equations, assuming the accretion flow is geometrically thin. The validity 
of the assumption that the disk is geometrically thin will be discussed later. As a first 
order approximation, we adopt isothermal structure in the vertical direction, which means 
that the sound velocity Cg = (p/p)^^^ is independent of z. Furthermore, we also assume 
radial and azimuthal velocities are independent of z. The vertical structure is obtained 
analytically by solving the remaining vertical component of Euler equation which represents 
hydrostatic balance in the vertical direction: 

^ = -P^ = -P^K z. (16) 

The last equation of eq. (|TB|) assumes that the disk is geometrically thin. The density 
distribution in the vertical direction is 

p{r,z) = p{r,0)exp(^-^^ . (17) 

Here, H is the vertical scale height defined by 

H = Cs/nK. (18) 

With the vertical structure given above, we can now integrate above basic equations in the 
vertical direction and rewrite them as follows. Mass conservation now gives 

M = 2nr^Vr, (19) 

where S (r) = p (r, z)dz = 2^n/2p (r, 0) H is the surface density. We have taken Vr to 

be positive for the inward flow so that the mass accretion rate M takes positive value. The 
height-integrated version of the radial component of eq. (^ is 

dVr , I dW 2\ W dlnflK 
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where W (r) = f^^oP (^5 = 2^71 /2p (r, 0) H is the height-integrated pressure. The last 
term of the right hand side of eq. (|20|) corresponds to the correction for the decrease of the 
radial component of gravitational force away from the equatorial plane (Matsumoto et al. 
1984). We can integrate the conservation of angular momentum (eq. p|) with respect to r 
as well as z and obtain following simple equation: 

M (/ - /in) = 27rr2aiy, (21) 

where /in is the specific angular momentum swallowed by the central black hole. The 
energy equations (eq. [|ll|, |jl2|) are vertically integrated together with the thermodynamic 



relations (eq. [jl^, [|l^) and now can be written as 

M{1- /in) -y- + 2TTrAie, (22) 




1) dr 2 (7 — 1) dr dr I dr 

7 + 1 dlnWe 37-1 dlni: d\nQ 



27rr (g-d - A.e) , (23) 



2(7-1) dr 2(7-1) dr dr J 

where Ajg = Xie {z)dz = y^ifAje (0) is the energy transfer rate from ions to electrons 
per unit surface area. Qrad radiative cooling rate per unit surface area which will be 

explicitly given in the next subsection. The energy equations can be written compactly as 
follows: 

Qa"dv. = Qtis - A.e, (24) 
Qadv,e = — Qrad- (25) 



2.3. Radiation Mechanism 

Following the work of Narayan & Yi (1995), we consider three processes for radiative 
cooling: bremsstrahlung, synchrotron radiation, and Comptonization of soft photons. In 
order to obtain the spectrum generated by the accretion flow and the global structure of 
the accretion disk, we need to solve global radiation transfer problem in the radial and the 
vertical direction which involves self-absorption and incoherent scatterings. We treat such 
complicated problem in a rather simplified way: 1) we assume locally plane parallel gas 
configuration at each radius and 2) we separate the Compton scattering process from the 
rest of radiation processes, i.e. emission and absorption. This will not make serious error 
because the gas is so tenuous that generated photons rarely experience multiple scatterings 
before they escape. Note that at such low frequencies as radio frequencies, the effect of 
free-free and synchrotron self-absorption is so large that the disk can not be considered 
optically thin, while the optical depth for scatterings is constant at all frequencies. We first 
estimate the spectrum of unscattered photons at given radius by solving radiative diffusion 
equation in the vertical direction. For the isothermal plane parallel gas atmosphere where 
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density configuration is given by eq. (0) the radiative diffusion equation can be solved and 
gives tlie energy flux F^, of the unscattered photons at given radius (see Appendix B): 



27r 
73 



5. 



1 - exp (-2V3r; 



(26) 



where r* = ^K,y(0)H is the optical depth for absorption of the accretion flow in the vertical 
direction with k,^{0) being the absorption coefficient on the equatorial plane. Assuming 
LTE, we can write = Xu/ {4:nB^) where Xu = Xi^.brems + Xi^.synch is the emissivity. Note that 
equation (p6|) includes the effects of free-free absorption and synchrotron self-absorption at 
low frequencies. 



2.3.1. Bremsstrahlung Emission 

In our model, electron temperature exceeds rest mass energy of an electron in some 
cases where electron-electron bremsstrahlung is important as well as electron-proton 
bremsstrahlung which dominates electron-electron bremsstrahlung in the classical 
temperature regime. Thus the cooling rate per unit volume for bremsstrahlung is 

Cr = + C- (27) 

Following Narayan & Yi (1995), we employ 

q~i = 1.25nlaTcafmeC^Fei {9e) , (28) 
where a/ is the fine-structure constant and is given by 

1/2 



1.781^,'''* for < 1, (29) 



^9 

Fei (^e) = ^ [In + 0.48) + 1.5] for 6, > 1, (30) 

ZTT 



and 



(31) 



X (l + l.We + el - 1.25ef/2) for Be < 1, 

q^^ = nlcrlafC^2Aee (In 1.1232^^ + 1.28) for Oe > 1. (32) 
With the cooling rate given above, we can write the emissivity Xi^,brems as 

Xu,hvems = q^.Gexp [^] , (33) 
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where G is the Gaunt factor which is given by (see Rybicki & Lightman 1979) 

(7^Af3^y^' for^<l (34) 



2.3.2. Synchrotron Emission 

For the relativistic temperatures for the electrons, and in the presence of magnetic field 
which is of the same order as equipartition magnetic field, synchrotron emission is also very 
important. Following Narayan & Yi (1995), the optically thin synchrotron emissivity by a 
relativistic Maxwellian distribution of electrons is 

WcH = 4.43 X 10-»^j^/' j , (36) 

where I'{x) is given by 



2.3.3. Compton Scattering 

We proceed to consider the effect of Compton scattering. We make use of the idea 
of energy enhancement factor which is derived by Dermer, Liang, & Canfield (1991) and 
modified in part by Esin et al. (1996). The energy enhancement factor rj is defined as the 
average energy boost of a photon. The prescription for rj is 

r] = cxp (s {A - 1)) [1 - P (j^ + 1, As)] + rj^^P {j^ + 1, s) , (38) 

where P is the incomplete gamma function and 



A = i + 49e + mi, 

(39) 



r7 — MZk 
'/max — , 

• _ lll(?;max) 

~ ln(A) • 
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is the optical depth for scattering: 

= 211^(7^11 X max f 1, — ) , (40) 

where Tes = (1 -\- nea^ / t^u)^^"^ is the effective optical depth (Rybicki & Lightman 
1979). Equation (^) gives the correct estimate for the optical depth for scattering in the 
presence of the absorption. Note that the simple treatment given above assumes that the 
cross-section is Thomson cross-section rather than exact Klein-Nishina cross-section. With 
the energy enhancement factor 7], the local radiative cooling rate is given by 

Q7ad = / duv (^) (41) 

We calculate the radiative cooling rate numerically and use it for the energy equation for 
electrons at each radius to calculate the global solution of the accretion flows. 



2.3.4- Calculation of the Spectrum 

We need more detailed specification for the calculation of the spectrum. Knowing the 
spectrum of unscattered photons and the probability that a photon will suffer scattering 
given times, we can calculate the Compton scattered spectrum if we know the information 
about how the original spectrum is modified after experiencing one scattering. Note that 
the average energy boost given in ( p6|) is the energy boost averaged over the Maxwellian 
distribution of electrons and can not be used directly for the spectrum calculations. Such 
problem was precisely discussed first by Jones (1968) and corrected afterwards by Coppi 
& Blandford (1990). We make use of the formula given by Coppi & Blandford (1990) and 
calculate the resulting spectrum. The remaining problem concerning Compton scattering 
is how to treat the spectrum of the saturated Comptonized photons. We assume that 
photons which are scattered more than times saturate and obey the Wien distribution 
oc exp {—hu/kTf,). 

Since large fraction of the emitted radiation is generated at the radius fairly close to 
the black holes in our model, we cannot ignore the effect of redshift due to the gravity and 
the gas motion. We include the gravitational redshift by simply taking the ratio of the 
energy of a photon when observed to its energy emitted at radius r to be y'l — r^/r. To 
treat redshift due to relativistic gas motion in a simple way, we concentrate on the face-on 
case where the optically thin assumption is most adequate. Thus we simply take the ratio 
of energy change for the redshift due to gas motion to be 1/J 1 — (w/c)^. 
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2.4. Numerical Procedure 

We solve numerically the set of equations given so far with the boundary conditions. 
The outer boundary conditions imposed are 

n = OmK, (42) 

Tgas=/i(- + -) =0.ir.ir, (43) 
Qrad = A.e (44) 

at Tout = lOOOOrg. Here Tvir is the virial temperature defined by 

Tvir = (7 - 1) (45) 

We did not set the angular velocity to be the Keplerian angular velocity itself at the 
outer boundary for simply technical reasons. We need to set somewhat sub-Keplerian 
disk at the outer boundary so that the viscous heating term take positive value with the 
simple viscosity prescription adopted in this paper [eq. (|10|) ]. We confirmed that the outer 
boundary condition have little effect on the structure of the inner advection-dominated 
flows. There have been many suggestions about how outer Keplerian disks are connected to 
the advection dominated disks. For instance, Honma (1996) took into account the effect of 
thermal conductivity and obtained the global solutions where the outer cool Keplerian disk 
is connected to the inner hot optically thin disk. However the mechanism of the transition 
is not yet clear. 

The free parameters in the set of equations are a, j3, Mbh, M, and li^. These five 
parameters are not independent because of the transonic nature of the equations, /jn can 
be determined uniquely so that the solution should satisfy the transonic condition when 
the remaining parameters are given. We have to adjust /in recursively to obtain smooth 
transonic solutions, since the location of the transonic point is not known until we obtain 
global transonic solutions. 



3. Results 

To discuss the general properties of our model, we set m = Mbh/Mq = 10^ and 
rh = M/Mc = 10~'^ where Mc = 32iTcrg/ Hcs and assign typical values for other parameters: 
a = 0.1, P = 0.5 (case of equipartition). To demonstrate how the heatings and coolings 
balance in the energy equations, we show Figure 1 the Q's for the ions (upper panel), and 
for the electrons (lower panel); (see eq. [^). The solid line in the upper panel shows the 
ratio of advective cooling of ions to the viscous heating, which is commonly denoted by /. 
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We have / = 1 for the purely advection dominated accretion flows and / = for the purely 

radiative cooling dominated accretion flows. Wc sec that / asymptotically approaches unity 
as the radius decreases. Wc use the region r < 2000rg to investigate the spectrum from the 
optically thin advection dominated accretion flows. We see from Figure 1 that the accretion 
flows becomes highly advection dominated in the region r < lOOrg where the heat transport 
from ions to electrons by the Coulomb coupling practically becomes zero. The interesting 
point is that it is electron advcctivc heating, rather than the heat supplied from the ions, 
that balances with the radiative cooling of electrons, which means the electron accretion 
flow becomes cooling flow near the central black hole where the electrons are cooled by 
consuming the stored internal energy to radiate before being heated up by the ions (cf. 
Nakamura et al. 1997). Many analyses concerning two-temperature accretion flows have 
adopted simplified energy equation for electrons: Ajg = Q^^d- results shows that above 
simplified energy equation is inadequate for the analysis of two-temperature advection 
dominated flows and one should adopt Q'^^ ^ = —Qra.d instead. 

Figure 2 shows the ion and electron temperature profiles. We see that the electron 
temperature is heated up to Tg ~ lO^'^K. Although we need to include the effect of 
electron-positron pair production and annihilation for such hot accretion flows, we make 
simple assumption that the pair density is very low (see discussion in Esin et al. 1996, 
Bjornsson et al. 1996, Kusunose & Mineshige 1996). 

Figure 3 shows the aspect ratio h/r of the accretion flow. We have used height- 
integrated equations for our calculations which neglect the higher order of h/r. We 
may conclude that the value oi h/r shown in Figure 3 is marginally safe for the height 
integration. However, of course, the vertical structure of the advection dominated accretion 
flows is an important issue and we will investigate it in future papers. 

Figure 4 shows the angular momentum and the various velocities as a function of the 
radius r. Here c* is deflned as 

/ (37-1) + 2 (7-1) a' W 
'^"V 7 + 1 S 

so that Vr — Cg at the critical point. We see from the Figure 4 that the radial dependence of 
the radial and the azimuthal velocities are different from that of the self-similar solutions, 
in which all the velocities are proportional to r~^/^, especially in the super-sonic region. We 
see that for the parameters given above, the azimuthal velocity v^p is fairly sub-Keplerian 
and even sub-sonic and of the same order as the radial velocity unlike the standard 
model, which is known to be a common feature for a = 0.1. (In the standard model, we 
have v^p ~ rQ^ = c^h/r ^ c* ^ Vr-) Note that we have 'sub-sonic' azimuthal velocity when 
the rotation is highly sub-Keplerian [v^p < rfl-K) and the disk is hot and geometrically thick 
{h/r ~ 1). For the reason that the radial velocity Vr is the same order as the azimuthal 
velocity v^, we call the accreting gas "accretion ffow" rather than to call "accretion disk" 
in this paper. 
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Figure 5 illustrates the luminosity distribution as a function of r. The contributions 
from respective radiation mechanisms (i.e. bremsstrahlung, synchrotron, Comptonization) 
are also shown. We see that the bremsstrahlung emission is important and the effect of the 
synchrtron emmision and the Comptonization is negligibly small in the outer region of the 
accretion flow (r > lOrg), while the synchrotron emission and the Comptonization dominate 
in the hot inner region (r < lOrg). The contribution of the Comptonization rapidly increases 
as the radius decreases because the amount of the synchrotron soft photons increases 
and not because the Compton |/-parameter y = TeskTe/nieC^ increases. Note that in the 
innermost region the electron temperature Tg is approximately constant and the surface 
density decreases and thereby y actually decreases. Note that almost all emission comes 
from the hot inner region (r < lOrg) and the fraction of the emission from the super-sonic 
region, which the self-similar solutions fail to describe accurately, is fairly large. As far as 
the bremsstrahlung is concerned, the contribution of the emission from the large radii is not 
negligible, which make the slope of the bremsstrahlung spectrum less steep, since what we 
observe is the superposition of the bremsstrahlung peaks from different radii. 

In Figure 6, we shows the spectrum generated by the optically thin accretion flows. 
Parameters are m — 10^, m — 10~^, a — 0.1, /3 = 0.5. S indicates the synchrotron peak 
which is composed of Rayleigh- Jeans slope and the optically thin synchrotron emission. CI 
and C2 indicate the once and twice Compton scattered photons, respectively. B indicates 
the bremsstrahlung emission plus photons suffering multiple Compton scattering. W 
indicates saturated Comptonized photons which form Wien tail. 

The upper panel of Figure 7 shows the surface density of the accretion flow. The 
dashed line corresponds to the case with the central black hole mass of Mbh = IOMq. 
We flnd that the structure of the accretion flow is nearly the same when the radius and 
the mass accretion rate are scaled with the Schwarzschild radius and the critical mass 
accretion rate, respectively. However we see in the lower panel of the Figure 7 that the 
electron temperature at the innermost region is varied slightly. This is due to the fact 
that the accretion flow is highly advection dominated. We remind the readers that the 
viscous heating is balanced almost entirely with the advective cooling in the ion energy 
equation and there is little coupling between the ions and the electrons. Thus the structure 
of the accretion flow is governed by ion energy equation while the electron energy equation 
including the radiative cooling is decoupled and determines the electron temperature. 
Although the height-integrated quantities are the same when the radius is scaled with 
the Schwarzschild radius, the amount of radiative cooling, which is the function of the 
density distribution rather than the height-integrated quantities, is different. Thus we have 
different electron temperatures for different black hole masses. We show the spectrum from 
the accretion flow for the IOM0 case in Figure 8. The shape of the spectrum is essentially 
the same as the lO^M© case, but the position of the synchrotron peak and the absolute 
luminosity differs considerably. 



- 15 - 



4. Application to Sagittarius A* 

Following up on the successful work of Narayan et al. (1995) which applied the 
advection dominated model to the Sgr A* (the central core of our Galaxy), we improve 
their model by fully solving the basic equations to explain the observed emission from 
radio frequencies to Gamma-rays. Figure 9 shows the model which explains observed 
radio and X-ray data quite well. The parameters assigned are given in the figure. The 
points and the short lines are the observational data, which assume interstellar absorption 
Nif = 6 X lO^^cm"^ and a distance d — 8.5kpc, both typical for the Galactic Genter, 
compiled by Narayan et al. (1995) (references therein). The various fines correspond to 
tfie different values of mass accretion rate wfiicfi varies by factor of 2. As Narayan et al. 
noticed, the position of the Rayleigh- Jeans slope is determined solely by the mass of the 
central black hole. Thus to give a good fit to the observed radio emission with our model, 
we have no choice but to fix the mass of the central black hole to be Mbh = 10^ Mq. 
Gonsidering the inaccuracy of the potential at the innermost region of the accretion fiow, 
this value is consistent with the value which is derived from the gas and stellar dynamics 
(Genzel & Townes 1987). We have M^. = 3.5 x lO~^M0/yr for the black hole of the mass 
-^BH = 10® Mq, hence the predicted mass accretion rate is ~ 2 — 5 x lO^^Mg/yr. 

Figure 9 also illustrates the M dependence of our model. The surface density sensitively 
depends on M, while the temperature is fairly insensitive to M. Thus the luminosity at all 
frequencies decreases when we reduce the mass accretion rate. We find that if we change 
the mass accretion rate by a factor of ~ 2, X-ray luminosity varies by the same factor 
while the synchrotron peak varies little. We suggest that the various X-ray data seemingly 
inconsistent with each other are purely due to the change of the mass accretion rate of 
the accretion flow. There have been an argument that the X-ray luminosity of Sgr A* is 
variable on the timescale of half a year, which is consistent with our suggestion. 

Figure 10 illustrates how the spectrum changes with different black hole masses. As we 
mentioned before, the structure of the accretion flow is the same with proper scalings, but 
the electron temperature is slightly different because of the difference of the emissivity. We 
see clearly that the position of the Rayleigh- Jeans slope is determined by the mass of the 
central black hole. 

We show in Figure 11 the /3-dependence of the structure of the accretion flow and the 
spectrum. We remind the readers that the parameter /? is defined as the ratio of the gas 
pressure to the total pressure. We have the stronger magnetic field for the smaller value 
of /3. Unlike the other parameter dependence, the ion temperature as well as the electron 
temperature decreases when we lower the value of f3. This is reasonable because for the 
lower value of /?, the stronger becomes the magnetic pressure and the gas pressure can be 
smaller to support the accretion fiow. This temperature change has a large effect on the 
bremsstrahlung spectrum but has little effect on the synchrotron peak, which is because the 
effects of the stronger (weaker) magnetic field and the lower (higher) temperature roughly 
cancel out. We set Tgas = 0.2Tvir at the outer boundary for P — 0.95 since we could not find 
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advection dominated solution for (3 — 0.95 with Tgas = O.lTyir. 

Narayan et al. (1995) did not investigate the parameter dependence of the viscosity 
parameter a and the mass accretion rate M independently, since they used simplified disk 
model, in which a and M cannot be chosen independently. Thus it is meaningful to study 
the a dependence of the spectrum. Wc show in Figure 12 the a-dependence of our model. 
When a is large, the angular momentum is extracted efficiently and the surface density 
decreases, which makes the bremsstrahlung emission weak. However, the increase of the 
electron temperature at the innermost region makes the synchrotron emission stronger. 
Wc find that the width of the synchrotron peak is the most sensitive to the value of a. 
The upper limits in the radio frequency band imposes strong restriction upon the value of 
a for the case of Sgr A*. For instance, we can not fit the entire spectrum with a — 0.1. 
We have to set a < 0.025 to explain entire spectrum of Sgr A* with our model. This fact 
is important because it is claimed that a ~ 0.1 in advection-dominatcd accretion flows 
on various grounds (see discussions in Narayan 1996). If we are to set a = 0.1, wc have 
to reduce the mass accretion rate by factor of 4 so as to give a good fit to the observed 
radio and IR spectrum (see Figure 12). In that case, we cannot explain the X-ray emission 
from the Galactic Center. Various X-ray observations have found the X-ray sources at the 
Galactic Genter but the angular resolution has been insufficient to identify a source with 
the radio source Sgr A*. There remains the possibility that X-rays do not come from Sgr 
A* at all (e.g. see Duschl et al. 1996). If that is the case, all the X-ray data are merely 
upper limits and the accretion flows with a = 0.1 do not conflict with the observations. 

There exist some data points in near IR and radio frequencies which cannot be 
accounted for. However, it is natural to think that there must be dust region or stellar 
contamination or non-thermal objects like jets in central region of our Galaxy overlapping 
the accreting black hole which in total we observe as a point source Sgr A*. In that sense, 
one should consider the observational data as upper limits. A point worth emphasizing here 
is that the entire spectrum of Sgr A* is basically explained with an optically thin accretion 
flow around a black hole of the mass Mbh = 10^ Mq. Moreover, we have solved globally the 
basic equations including the gas dynamics and the radiation processes consistently rather 
than to consider more primitive models like isothermal gas complex of certain size. 



5. Summciry and Discussion 

In this paper, we have calculated the global structure of optically thin advection 
dominated accretion flows in the context of two-temperature plasma, adopting the full set 
of basic equations including the energy equation for the electrons. We have also calculated 
the spectra emitted by the optically thin accretion flows which we calculated. Wc have 
made the calculation of the spectra and that of the structure of the accretion fiows to be 
completely consistent by calculating the radiative cooling rate at each radius by numerically 
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integrating the whole spectrum emitted at the radius. 

As a result of the advection domination for the ions, the heat transport from the 
ions to the electrons becomes practically zero and the radiative cooling balances with the 
advective heating of the electrons. This means that the electron cools itself by releasing 
the stored internal energy as a radiation. Hence the energy equation for the electrons play 
an important role for the calculation of the spectra, where the temperature profile of the 
electron is the important factor. 

An interesting feature of the advection dominated flow, which is known already, is 
that the azimuthal velocity becomes highly sub-Keplerian and of the same order as the 
radial velocity and the sound velocity. The point worth noting is that in the innermost hot 
luminous region, the divergence of the velocities from those in the self-similar solution is 
fairly large. 

The accreting gas becomes very hot. The electron temperature even exceeds the 
rest mass energy of an electron. We have not taken into account the effect of the pair 
production and the annihilation, which is an important issue. For such hot accretion flows, 
the synchrotron emission and the Compton scattering are very important. 

The spectrum is composed by 1) the synchrotron peak which comprises optically thin 
synchrotron emission and the self-absorbed Rayleigh- Jeans slope and 2) the unsaturated 
Comptonized photons which forms some bumps and 3) the bremsstrahlung emission and 
4)the saturated Comptonized photons. The dependence of the each component on the 
model parameters is complex. Among them, the position of the Rayleigh- Jeans slope 
is almost solely determined by the mass of the central black hole. When we make the 
magnetic fleld stronger, the temperature of the entire flow decreases, which has signiflcant 
effect on the Comptonization and the bremsstrahlung emission, but has little effect on 
the synchrotron emission. When we make the viscosity smaller, the surface density 
increases and the bremsstrahlung emission increases, but the synchrotron emission and the 
Comptonization decreases. The simplest relation is the dependence on the mass accretion 
rate. When we reduce the mass accretion rate, the entire emission is reduced. However the 
bremsstrahlung emission is much more sensitive to the change of the mass accretion rate 
than the synchrotron emission. 

We find that the spectrum of Sgr A* is explained by the optically thin advection 
dominated accretion flow around a black hole of the mass Mbh = IO^Mq. Narayan et al. 
(1995) also calculated the spectrum of Sgr A* using an optically thin advection dominated 
accretion flow model. Their best flt parameters are different from ours. For instance the 
mass of the central black hole is Mbh = 7.0 x 10^ Mq according to their model. The different 
points in our model are 1) full global treatment of the basic equations and 2) inclusion of 
the electron energy equation and 3) calculation of the innermost region where the flow is 
supersonic and the effects of the redshift are important. 2) and 3) have very important 
effect on the emitted spectrum, which is not considered in Narayan et al. (1995). We 
conclude that the X-ray data obtained by various satellite observations are explained by the 
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variation of the mass accretion rate by a factor of ~ 2, if we allow a to have small value. 
If we set a ~ 0.1, which is considered to be a standard value for the advection dominated 
accretion flows, it is not likely that the X-rays come from Sgr A*. 

We have computed the model using height-integrated equations with fixed structure in 
the vertical direction. We have also adopted simplified form of equations for the radiation 
field. Our immediate goal is to solve the basic equations including equations for the 
radiation field in at least two dimensional space. However, our successful result presented 
in this paper tells us that the basic idea is correct. The full treatment of Schwarzschild or 
Kerr metric is also an important issue. 

We thank Professor Shoji Kato for useful discussions, and Ann Esin and Jun Fukue 
for useful comments. We are very grateful to R. Narayan for providing us with a code to 
calculate inverse-Compton spectra. We thank the referee J. P. Lasota for many valuable 
comments and suggestions. 



In this paper, we have adopted the height-integrated equations to obtain the structure 
of the optically thin advection dominated accretion fiows. We have seen that the 
temperature of the accretion fiows is very high and the aspect ratio /i/r is ~ 0.5. Thus 
the height-integration may not be an excellent approximation since the operation of the 
height-integration omits the higher orders of h/r. Thus wc consider it useful to describe in 
detail the operation of the height-integration to discuss the reliability of our calculations. 

The potential we adopted is spherically symmetric, while we adopt the cylindrical 
coordinate system. Thus we need to expand it around the equatorial plane {z = 0) in 
order to integrate the basic equations in the vertical direction. In the case of geometrically 
thin accretion flows, we can neglect the terms containing the higher orders oi z/r and the 
operation of the height-integration becomes very simple. 

The potential is expanded as 



A. Height Integrations of the Basic Equations 



V'(r,0) 1 + 



1 d\n%l) 

2 dhir 



z=0 




(Al) 



Thus we approximate as 



(A2) 



where Jlfe is defined as 




(A3) 
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Using Equation (A.l), the derivatives of are also expanded as follows: 

/ -v4' 



dxjj dip 
dr dr 



2=0 



amr r"^ 



dtp dip 
dz dr 

We approximate the derivatives of ip as follows: 



z 








r 

z=0 





dr \ dlnr r'^ 



^ = ah- 

dz " 



(A4) 
(A5) 

(A6) 
(AT) 



We have used eq. ( |A^ ) to determine the vertical structure of the accretion flows [see eq. 

mi 

We then integrate the basic equations in the vertical direction. The integrations of the 
continuity equation (eq. 0) and the azimuthal component of the Euler equation (eq. 0) 
are straightforward. In the process of the integration of the radial component of the Euler 
equation, we encounter the following integration: 



J-oo dr " 



W- 



dlnQk 



dr 



(At 



The last term in the right hand side of eq. ( [A8| ) is sometimes omitted in the height-integrated 
equations which appear in the papers on the accretion disks, which is not proper, since it 
can be the same order as the other terms in the Euler equation. 



The left hand sides of the energy equations (eq. |jTT|], eq. 0]) together with the 
thermodynamic relations (eq. |]T^, eq. |TB|) are transformed to give 



/ ds ds^ 
\ dr dz , 



We can integrate Equation 



7 



Id, , d , , 

-— {rVrP) + TT- [VzP) 

r dr dz 



7-1 
with the integration 



dp 
v^—dz 
dz 



-WVr 



dlnH 



dr 



dp dp 



(AlO) 



to obtain the height-integrated version of the Energy equations (eq. ||2^, p3|). 
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B. Calculation of The Flux of Unscattered Photons 

Consider the isothermal plane parallel atmosphere with the density distribution being 



p{z)= p (0) exp 



(Bi: 



Assuming the Eddington Approximation which is valid for isotropic radiation fields (and 
even for slightly nonisotropic fields, see Rybicki & Lightman 1979), the radiation field in 
the vertical direction is described by the radiative diffusion equation: 



3 



(B2) 



where Ti, is the optical depth from the surface of the accretion fiow. There is no well-defined 
surface of the accretion fiow since the density tends to zero when we increase z but will 
never equals to zero [see eq. (|Bl|)]. However, the optical depth of the accretion fiow in 
the vertical direction is finite and hence we can define the surface when the vertical height 
is measured with the optical depth. Note that t = t* = ^k,^{0)H at the equatorial plane 



and r = 2r* at the other surface. We solve eq. 



with boundary conditions. We take 



1 dJ,, _ J 

lt = (r. = r;). 



0), 



(B3) 



The boundary condition at the surface assumes there is no irradiation onto the surface of 
the accretion fiow and is derived by adopting two-stream approximation (see Rybicki & 
Lightman 1979). The solution for eq. ( p2D which satisfies the boundary conditions [eq. 
is 

p—^Ti, 



J. 



B,A1 



(B4) 



The energy fiux on the surface of the accretion fiow is given by 



47r dJ^ 



3 dr,. 



7! ^ 



1 - exp (-2V3r; 



(B5) 



Note that 



FAO) = f^B, (r:»l), 

F,{0) = ^Xu{z = 0)H (r;«l) 



(B6) 
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Fig. 1. — The amount of the heating and the coohng rates for the ions (upper panel) and for 
the electrons (lower panel) in the units of 2T^r'^ / {MWi/T) and 27rr^/(MWe/S), respectively. 
The ratio of advective cooling of ions to the viscous heating / is also shown in the upper 
panel. The accretion flow becomes advection dominated (/ = 1) as the radius decreases and 
the energy transport from the ions to the electrons becomes practically zero. The radiative 
cooling is balanced with the electron advective heating in the innermost region. 

Fig. 2. — The temperature of the ions and the electrons as a function of r. The electrons 
attain the temperature of 10^°K. 

Fig. 3. — The aspect ratio h/r of the accretion flow. 

Fig. 4. — The angular momentum (upper panel) and the various velocities (lower panel) of 
the accretion flow. The azimuthal velocity is highly sub-Keplerian. In the hot luminous 
region near the central black hole, the divergence of the velocities from those in the self- 
similar solution is fairly large. 

Fig. 5. — The luminosity distribution as a function of r. The contributions from respective 
radiation mechanisms (i.e. bremsstrahlung, synchrotron, Comptonization) are also shown. 
The synchrotron emission and the Compton scattering dominate in the innermost region. 
The distribution of the bremsstrahlung emission is relatively flat. 

Fig. 6. — The spectrum generated by the optically thin accretion flow around the central 
black hole of the mass Mbh = lO^M©. S indicates the synchrotron peak which comprises 
Rayleigh- Jeans slope and the optically thin synchrotron emission. CI and C2 indicate the 
once and twice Compton scattered photons, respectively and B indicates the bremsstrahlung 
emission plus photons suffering multiple Compton scattering. W indicates saturated 
Comptonized photons. 
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Fig. 7. — The surface density (upper panel) and the temperature (lower panel) of the 
accretion flow. The dashed hne corresponds to the case where Mbh = WMq. We find 
that the structure of the accretion flow is nearly the same when the radius and the mass 
accretion rate are scaled with the Schwarzschild radius and the critical mass accretion rate, 
respectively, while the electron temperature at the innermost region is varied slightly. 

Fig. 8. — The spectrum generated by the optically thin accretion flow around the central 
black hole of the mass Mbh = lOM©. The letters indicate the same meaning as Fig. 6. The 
shape of the spectrum is essentially the same. However, the position of the synchrotron peak 
and the total luminosity differ significantly. 

Fig. 9. — The spectrum of Sgr A* (upper panel). The lines correspond to the spectrum 
calculated with our model presented in this paper. The M-dependence of the spectrum and 
the temperature and the surface density (lower panel) are illustrated. The spectra behaves 
in the simplest way. When M is reduced, the surface density and the entire emission are 
also reduced. The bremsstrahlung emission is more sensitive to the M change than the 
synchrotron emission. 

Fig. 10. — The MeH-dependence of the spectrum (upper panel) the temperature and the 
surface density (lower panel). The position of the Rayleigh- Jeans slope determines the mass 
of the central black hole. The surface density is not varied when Mbh is changed. 

Fig. 11. — The /3-dependence of the spectrum (upper panel) the temperature and the surface 
density (lower panel). When (3 is lowered, the temperature decreases and the bremsstrahlung 
emission is weakened. However, the synchrotron emission is fairly insensitive to the change 
of /3. 

Fig. 12. — The a-dependence of the spectrum (upper panel) the temperature and the 
surface density (lower panel). When a is lowered, the surface density increases and so the 



-25- 



bremsstrahlung emission. However, the synchrotron emission weakened due to the decrease 
of the electron temperature in the innermost region. 
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